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Abstract 

The heat transport in highly porous fiber structures is investigated. The fibers are supposed to 
be thin, but long, so that the number of the inter-fiber connections along each fiber is large. We 
show that the effective conductivity of such structures can be found from the correlation length 
of the two-point correlation function of the local conductivities. Estimation of the parameters, 
determining the conductivity, from the 2D images of the structures is analyzed. 
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I. INTRODUCTION 


The composite materials comprising long libers of high thermal conductivity immersed in 
some amorphous matrix of low thermal conductivity are now quite popular in applications. 
Such materials can be used as the thermal barrier coatings. The long fibers are necessary 
to ensure the mechanical strength of the coating. Long high-TC fibers are also useful to 
enable anisotropic properties of the coating. Ablative thermally protective systems is an 
example of such structures. These systems typically consist of low density but rigid carbon 
fiber material infiltrated with the phenolic resin. Such structures have a very low density 
(~ 0.27g/cm 3 ; this means that fibers occupy approx. 12% of the bulk of the specimen), 1 the 
typical diameter of the carbon fibers is 14-16 //.m and their length exceeds 1600 /im. Using 
special manufacturing methods enables the parallel orientation of fibers in order to promote 
minimum through-the-thickness thermal conductivity, which is one of the critical properties. 

In spite of many detailed studies, potential of highly-porous fiber materials seems to be 
not realized. Further progress depends on the possibility of accurate estimation and mod- 
eling of their properties. There exist some sophisticated approaches and software packages 
for modeling of these protection systems, for example the program package FIAT (Fully Im- 
plicit Ablation and Thermal response code) developed by NASA for one-dimensional (ID) 
simulation of kinetically controlled pyrolysis, in-depth conduction, blowing due to pyroly- 
sis gases, surface recession as the function of time using ID model of ablation of porous 
material. 2-4 The program package TITAN, also developed by NASA, allows 2D modeling. 5 
Several other software packages for modeling of fiber materials properties are available: TI- 
TAN/GIANTS (integrated procedure), MARC (commercial FEM code) 6 , and data base 
TPSX (http:/ /tpsx.arc. nasa.gov/). 7 

Nevertheless, direct numerical calculation of the transport properties in 3D fiber compos- 
ites is very difficult (see, e.g., Ref. 8). Indeed, it is necessary to consider the 3D transport 
equation in a fine-structure mesh. The mesh spacing should be much smaller than the fiber 
thickness, while the total size of the mesh network should exceed the fiber length. More- 
over, the composite structure is very complicated and its microscopic description requires 
a lot of parameters. 9 The question arises whether all of these parameters are important for 
the thermal conductivity (TC) estimation or whether it can be estimated from a few other 
parameters related to them. We have shown that the two-point correlation function (CF) of 
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the local conductivities is the only parameter, which is necessary for TC estimation in the 
non-homogeneous materials of low porosity. 10,11 It is not clear, however, whether the CF is 
sufficient for the modeling of thermal conductivity of highly-porous fiber materials (porosity 
~ 80 - 90%) or whether the higher order correlation functions should also be employed. 

An important property of the fiber structures is their anisotropy. Anisotropy can be en- 
abled artificially, however, a typical fiber structure is always anisotropic due to the nematic 
ordering. Effect of anisotropy on the electrical conductivity in such structures has been 
studied 12 and exact scaling laws have been obtained. 13 The dependence of thermal conduc- 
tivity of fiber composites on the fiber orientation was observed, e.g., in oriented fibrous 
carbon insulation. 2 The room-temperature in-plane thermal conductivity of the carbon in- 
sulation at 1575 K was 0.05 W/mK normal to the fibers plane and 0.14 W/mK in the 
fibers plane. It was also reported, that the thermal conductivity of the fiber materials of the 
phenolic resin infiltration in the direction of fibers alignment is ~ 2.4 higher than that across 
the thickness. 14 Also filled honeycombs, multilayers, tile-like structures are manufactured to 
provide anisotropic properties of the structure. 15 

Qualitative estimation of conductivity of the fiber structures has been done in Ref. 16. 
It was shown that the conductivity can be expressed via porosity of the structure, thickness 
of the fibers, and mean distance between fiber connections. It is not clear, however, how to 
estimate these parameters from the two-point correlation function, which is not simple since 
CF cannot be approximated with a single exponent or Gaussian function. This problem 
will be considered in this paper. We show that total porosity and the correlation lengths 
along each (x, y, and z) direction are the only parameters necessary for estimation of the 
conductivity of the fiber structures. 

We start with a simple model of TC, which is a generalization of that proposed in Ref. 16 
for an anisotropic media. We suppose the fibers are long enough, so that each of them 
has many crossing points along its length. This means also that the number of fibers per 
unit volume essentially exceeds the percolation threshold. This is the main point which 
distinguish our model from the previous studies where the conductivity near the percolation 
threshold has been investigated. 12 We have found that TC can be expressed via porosity 
of the structure and mean distance between the fibers in each direction. In the following 
two sections we show that all these parameters can be estimated from the CF of the local 
conductivities. 
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In Sec. Ill we consider the 2D models of the structures consisting of straight square 
fibers of equal thickness and length. By changing the length of the fibers, we can measure 
all parameters of the structure and compare them with the CF calculated using the image 
of the structure. We relate the parameters of the structure with the correlation lengths of 
CF. We show that CFs of the structures with the same fiber thickness and main distance 
between the fibers are close. In Sec. IV we generalize this result to the 3D case. 


II. THERMAL CONDUCTIVITY OF FIBER STRUCTURES 

A. Thermal conductivity of the Fiberform 

Fiberform is the structure composed exclusively of fibers without the resin infiltration. 
The heat propagation through the fiber network is the only mechanism of the thermal 
transport in such structure. By assuming the Fiberform as an array of cylindrically shaped 
fibers, each fiber can be geometrically characterized by its length l and diameter d. If n is the 
fiber concentration, i.e. the number of fibers in the unit cube, then these three parameters: 
l, d, and n are the characteristic of the Fiberform. We introduce also the Fiberform density 
as Vf = Tind 2 l/A (V/ = ndl in the 2D case). All fibers are intersecting in a 2D model, if they 
are long enough and close. This is not the case in the 3D model, as not coplanar fibers can 
be skewed. 

We suppose that the fibers are thin (d <C /) and long, so that the mean distance between 
the neighboring connecting points A is small compared to the fiber length, (A <C l). This 
also means that the concentration of the fibers is large enough, so that we are well above 
the percolation threshold. 

Suppose an external thermal gradient applied along the 2-axis, and assume the cylindric 
symmetry of the Fiberform with regards to this axis, i.e., we consider the heat flow along 
the 2-direction (Fig. 1). Let us choose two points in the fibers at the plane 2 = const. 
The temperature, T, at these points can be different but its average values are equal. In- 
deed, the condition A l means the existence of percolation between the points, while 
( T{ri )) 7^ (T(rfc)) ensures the heat flux between them. The latter is forbidden by the 
cylindric symmetry, therefore (T(rj)) = ( T(rk ))• 

Let us find the heat flux, Iij, through a fiber connecting two points r* and Vj at two close 
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FIG. 1: Schematic of heat transport in a fiber structure. 2 is the direction of heat flux, 
planes z = 0 and z = h, respectively (see Fig. 1): 


li i 


FW - T(r, 


Here Kf is the fiber thermal conductivity, / y - = |r* — r 7 1 = /i/| costal, and is the angle 
between the fiber and the 2 -axis. The total flux through the plane z = 0 is equal to sum of 


all fluxes over all the fibers crossing the plane 




[{T(r,)> - <T(r,)}]. 


This allows us to determine the effective thermal conductivity /r, ( . ir as: 


Keff “ 4 S[{T(n)) - (T(r,))] ( cos0 ^T(ri) T (u)])- 

Here N = AS/ kD\ is the total number of fibers crossing the plane z = 0, S is the side area 
of the specimen in this plane, and D± is an average distance between the fibers crossing the 
plane z = 0. Supposing that (cos [T(r, : ) — T(rj)]) = (cos djj}([T(rj) — T(rj)]), one finds 

d 2 .. _ .. D\\ 

Keff _ D 2 01 — 2D|| + ^ 

where D\\ is an average distance between the fibers along the 2 -direction, so that Vj = 
d 2 (2D± + D\\)/ 

In a similar way we can consider the TC across the 2 -axis. The results are 


K \\ jj 2 lx 'f ’ Kj - 
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or equivalently, 


D I, 


Kll = 


-VfKf, K _ L - 


D\ 


( 2 ) 


2D|| + D. J •' 2D\\ + D 

where k\\ and are the components of the effective TC along and across the z-axis, respec 
tively. Thus, 


K± + 2/C|| — VfKf, (3) 

/' I //' = T>J_/ D|* 

The former equality represents the trace of the conductivity tensor, which is invariant with 
respect to the axis used. The latter equality is in agreement with Ref. 13. In particular, for 
an isotropic structure D\\ = D± = D and the effective TC becomes 

_i T/ 

K eff — jj2 K f ~ 3 VfKf. 

The factor 1/3 appearing in the last equation indicates that the rule of mixtures cannot 
be applied to the fiber structures. To understand the reason, let us consider the TC along the 
side of a simple cubic fiber network. We assume the thin fibers, so that the heat propagation 
exists only along the fibers, but does not across them. Therefore, only the fibers oriented 
along the applied temperature gradient participate in the TC. Effect of the other 2/3 fibers 
can be estimated, if we replace the cube side D by D — d. Then = d 2 / (D — d)‘ 2 ~ 
(1 + 2 d/ D)d 2 / D 2 . The factor 1/3 disappears when d ~ D. 


B. Thermal conductivity of ATPS 

ATPS is the material comprising the fibers immersed in some matrix of small TC, km 
K f. The effect of the matrix can be two-fold. First, the direct heat flux across the matrix. 
Second, the heat transport through fibers and narrow bridges between two skewed fibers. 
The second effect can be important, if thermal resistivity of such bridges is of about the 
resistivity of the fibers between two connection points. Let us first consider the second effect. 

Let 5 be the length of such a bridge. The bridge is important if its resistivity is less 
or about the resistivity of the appropriate fiber part, i.e. 8/km 1$ A*//'/- Here A* is the 
distance between two bridges which size is less than 5. It can be estimated similar to that 
of Eq. (9): 

3 

2 n(d + 5/2)1 
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FIG. 2: Two dimensional models of fiber structures and corresponding correlation functions, W, 
taken along the horizontal axis, x. (a) Short-fiber models with aspect ratio, (l/d) = 2,6 and 10, 
Vf = 0.01, 0.04, and 0.07, respectively, (b) Long-fiber models with aspect ratio, (l/d) = 10,20 and 
36, Vf = 0.13, 0.17, and 0.21, respectively. All fibers have the same positions, orientations and 
thickness in each model. A thickness value of d = 0.5 unit was used in the computations. 


Considering the cube, A* x A* x A*, we can write its resistivity as 


R 


A* S \ D 2 _ 1 

+ W 
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FIG. 3: Long fibers. Small circle intersects the fiber two times (bold line). Large circle has many 
points of intersection (dashed line). 


from which one finds: 


^eff 


\*d 2 /D 2 


d 2 /D 2 


( 4 ) 


A */ Kf + 5/ Km 1/kj + 5 / X*km 

The above relation coincides with a similar estimation for K e rr in Ref. 16, if we suppose, 
5/X* = d 2 / A 2 or 5 — Vfd. Equation (4) turns into Eq. (1), if the last term in its denominator 
is omitted. This is possible when 5/X* <C k,m/k / or Vf5 2 /d 2 <C k«/k/. 

Another way to take the bridges into account consists in assuming each fiber surrounded 
by some layer, so that the thickness of each fiber is d 0 > d and its effective TC is, k = 
[i Kfd 2 + Km(5q — d 2 )]/^. Substituting the latter into Eq. (1) with d —> do yields: 

K f d 2 + k m {51 - d 2 ) 


/’Vll 


D 2 


This means that the fiber network and the matrix can be considered as two resistors con- 
nected in parallel. Therefore, the TC of the structure can be estimated as: 17 


1 4 

K eff = gCf'Af + (1 - 


( 5 ) 


III. CORRELATION FUNCTION OF THE 2D FIBER STRUCTURES 


In order to estimate the parameters D\\ and D± required for the TC modeling, we consider 
the CF of fiber structures. We introduce CF as 


W(r,r') = 


(Mr) - (r/(r))][//(r / ) - (rj(r J ))]) 
4k/(l — Vf) 


( 6 ) 
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where (• • • ) denotes an ensemble average and the characteristic function rj(r) is given by 


r/(r) 


1, if r is inside a fiber 
— 1, otherwise. 


For the statistically homogeneous media, Eq. (6) depends only on the difference between 
two arbitrary points r and r', [W (r, r') = W(r — r')]; it is equal to unity at the coordinates 
origin, fF(0) = 1 and vanishes at the infinity. If the media is also isotropic, then bF(ryr') = 
W{\r — r'|). The aim of this section is to calculate the CF using the structure image and 
estimate the parameters of the structure from the appropriate correlation lengths. We will 
see that the CF is determined by the fiber thickness d and the distance between the inter- 
crossing points of the fibers, if their density is large enough. Otherwise, the fiber length 
becomes the second correlation length. In evaluating CF from digital images we employed 
the procedure outlined in Ref. 18. 

Figure 2 displays the 2D fiber models and the corresponding CF. All the structures 
are composed of similar fibers; they have the same positions and orientations, and can be 
distinguished by the fiber length only. Fibers in Figure 2(a) are short, so that they do not 
intersect. The corresponding CF have two correlation lengths: one of them is relevant to 
the fiber thickness (r/d = 1) and the other one to the fiber length. The CF of different 
structures of this figure are well distinguished. 

On the contrary, the CF in Fig. 2(b) are very close. This is because the second correlation 
length for the long fibers is the mean distance between the inter-crossing points, but not 
the fiber length. Direct calculation of this distance (from the image using, e.g., Photoshop 
softwear) yields about 9 d for all the structures in Fig. 2(b). 

In order to understand the difference between the CF behavior in Fig. 2(a) and Fig. 2(b), 
let us consider the two-point probability function 19,20 


S = (rjo(r )//o(r')) , 


( 7 ) 


Vo (r) = 


where 

r 

1, if r is inside a fiber 
0, otherwise. 

Then S = Vf, if | r — r'\ d, i.e. if both r and r' are either inside or outside the 
liber. In order to estimate S for a longer distance, \r — r'\ ^$> d, we have to multiply the 
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FIG. 4: Cross-sections images and corresponding correlation functions of 3D models of fiber struc- 
tures. (a) Fibers randomly positioned and oriented in the space, (b) Cross section image of (a) 
at £ = 15. (c) Correlation functions W(x, 0) and W(0,y) of (b), taken along the x and y axis 
respectively, (d) Fibers randomly positioned with preferred orientation along the z-axis. (e) Cross 
section image of (d) at x = 10. (f) Correlation functions of (e) taken along the y and £ axis. A 
total number of N = 2000 cylindrical shaped fibers with aspect ratio (l/d) = 40 and Vf = 0.21, 
were considered in the simulations. 

probability that the point r belongs to the fiber [P(r E fiber) = Vf] by the probability for 
the point r' belongs to the fiber [P{r' E fiber)] provided by the first point r is in the fiber. 
The latter probability can be easily estimated, if we suppose isotropy of the fiber structure 
and draw a circle around r, as shown in Figure 3. Then P(r f E fiber) can be estimated 
as the ratio of the fiber thickness inside the circle to the circle length. For the short fibers 
[Fig. 2(a)], P(r f E fiber) = d/7i\r — r f |, if d | r — r'\ Z, and P(r f E fiber) = V/, 
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if | r — r'\ l ; in the latter case both probabilities are independent. For the long fibers 

[Fig. 2(b)], P(r' G hber) = d/ir\r — r'|, if d <C |r — r'\ l c , where l c is the mean distance 
between the crossing points. Note that P(r' G fiber) has an increase, when | r — r'\ > l c . 
Such increase can be found in the CF of Fig. 2(b) at | r — r'\ ~ lOd. For | r — r'\ / c , 

S = Vf as P(r G fiber) and P{r' G fiber) are independent at the large distances. From the 
above, Eq. (7) can be rewritten as 


Vf, 

jr — r'\ 

d, 

d 

d <C \r 

- r'\ 

n\r—r'\ 5 


V 2 
v f ’ 

\r — r'\ 

| > lc- 


5 = 

Since r/(r) = 2// 0 (r) — 1, it follows from Eqs. (6) and (7) that 

W(r - r') = 

thus 

W(r - r') = 


S{r - r') - Vf 

n(i-n) ’ 

| r — r'\ d, 


n\r—r'\(l — Vf) 

o, 


d < \r- r'\ < l c , 


( 8 ) 


|r - r'| > l c . 

A similar behavior of CF was found in Fig. 2(b). Note that d/ \r — r'\ A> V/, if |r — r'\ <C l c - 
Apparently, it is the connectivity (A), i.e. the mean distance between the crossing points, 
but not the mean length of the fibers that determines the transport properties in the fiber 
structure. We have shown that the connectivity coincides with the correlation length of the 
CF in the 2D case, A ~ l c . Therefore it can be estimated from Eq. (6). It is the connectivity 
which determines the conduction pathways. They are roughly the same in the Fig. 2(b) thus, 
all these structures have the same TC. Elongation of the fibers in this case means increase of 
the dead ends, but it does not change the conductivity. Increase of the TC can be expected 
when the fibers becomes so long that the number of the crossing points increases. This 
means decrease the connectivity value and, therefore, increase of the conductivity. 


IV. CORRELATION FUNCTION OF THE 3D FIBER STRUCTURES 

The main difference between the 2D and 3D models concerns the intersection of the 
fibers. Indeed, the fibers in the 2D model intersect, if they are long enough and close. In 
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general, the fibers in the 3D structure can be skewed; this happens with the fibers which 
are not coplanar. If so, then the mean distance between the crossing points should exceed 
the distance between the fibers. The question arises what is the second correlation length 
discussed in the previous section? Is it relevant to the connectivity or to the mean distance 
between the fibers? 


Let us first estimate the connectivity A. Consider a large parallelepiped, L x L x A 
(. L — > oo). The number of fibers within it is N — nL 2 X, and their total projection on the 
L x L plane is, (2/3 dl)N. This projection should cover the entire plane, if A is the distance 
between the crossing points. Thus, 

a = 2L 

2ndl 

The mean distance between the fibers D , can be estimated as D 2 nl = 1 , or 


Then for the ratio A/D, we obtain 



a _ 3 nr i 

where Vf = nnd 2 l / 4 is the density of the 3D fiber structure. This value is not so large for 
Vf ~ 0.1. 

It is interesting to compare A with the fiber length l : 


A _ 3 

7 “ 2 mH 2 ' 


From this relation it is evident that A l, if ndl 2 ~ 1, i.e., the fiber concentration 
n c = (dl 2 )~ 1 determines the percolation threshold. Indeed, at n = n c we have only one fiber 
in the l x l x d box. The fibers are not connected, if n <C n c , but connected in the opposite 
limit. The mean number of connections per fiber is, l/D ~ In this paper we assume 

this value to be large. 

For a rough estimation of the correlation function in 3D we have to replace the circle in 
Fig. 3 with a sphere, so that the conditional probability P(r' G fiber) provided by the first 
point r is in the fiber is now, P(r' G fiber) = d 2 / 8 | r — r ' | 2 . Therefore Eq. (8) becomes: 


W(r - r') 


1 , 

8 |r— r'| 2 (l— Vf) 

0, 


| r — r' | d, 

~T—Vf ’ d < l r - r '\ < 

I r — r ' I S> D. 


( 10 ) 
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(b) 



(c) 



FIG. 5: Cross-section of a 3D Fiberform (a), its correlation function, IF (b), and Fourier transform 
of IF in a; and y directions. [Image (a) provided to NASA courtesy of Xradia, Inc.] 

For the correlation length we obtain thus, l c ~ D. 

Just with the purpose to illustrate how the 3D case can be treated, we have simulated 3D 
fiber structures and extracted the cross-section images from them (Fig. 4). The algorithm 
employed to produce such structures allows free intercrossing of the fibers, so that nematic 
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ordering does not occur. It relies on the basic ideas exposed by Berryman, 18 but it has 
been readapted to the cylindrical geometry in this work. Figures 4(b) and 4(c) displays one 
of such cross-sections (z = 15) and the corresponding W(x, 0), TV(0, y) of a fully random 
fiber structure. As it can be appreciated, these functions are very close because both x 
and y direction are equal. The CF allows to estimate the mean distance among fibers as 
D = 2d. This does not occur when the fibers have a preferred orientation [Fig. 4(d)] and we 
examinate a cross-section containing the preferred axis [Fig. 4(e)]. The CF measured along 
the preferred axis and along the direction perpendicular to it, display different correlation 
lengths as shown in Figure 4(f). 

It is interesting to consider the Fourier transform of the correlation function. For the 
simple CF given in Eq. (10) and d <C D, it has the form 


(1 - V,)V 


(sin kD — kD cos kD ) . 


In Eq. 11 we omit the small values, which oscillate as cos kd. Thus, we can estimate the D 
value analyzing the oscillations of Fourier transform of the CF. This is important for the 
correlation functions measured from the structure images. 

Figure 5 illustrates the cross-section of a 3D Fiberform material, its correlation function 
and the Fourier transform. It is not easy to estimate the correlation length directly from 
the CF, because of the noise [Fig. 5(b)], However, it can be estimated from the Fourier 
transform: the period of the shortest oscillations is equal to 2i t/D. For an anisotropic 
structure, as that shown in Fig. 5(a), this criteria allows estimation of the ratio D±/D\\, 
which determines the anisotropy of the thermal conductivity along and across the structure 


axis. 

Let us estimate the thermal conductivity of the Fiberform presented in Fig. 5. Volume 
fraction of the fibers can be calculated directly from the digital image of Fig. 5(a), we found 
it equal to Vf = 0.1. Diameter of the fibers is the first correlation length of CF; we found 
d ~ 15 pm from Fig. 5(b). If we assume for the fiber length l = 1600pm, then for the 
fiber concentration, n = AVf /ird 2 l : we found n ~ 3.5 • 10 5 cm -3 ; this value much exceeds the 
percolation threshold n c = 1/dl 2 ^ 2.6 • 10 4 cm -3 . Rough estimation of the mean distance 
between the fiber connections is, A = 3/2 ndl ~ 180 pm; so that there are about 1 / A ~ 10 
connections for each fiber. Thus, all assumptions previously made to estimate the heat 
transport are satisfied. The distances between two adjacent peaks in both, bold and dashed, 
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curves in Fig. 5(c) are equal to 0.045 jum 1 . This means that Dg = D± ~ 140 /mi. This 
estimate can be observed in Fig. 5(b) as well. Then Eq. (3) yields k\\ = n± ~ 0.03 Kf. 


V. DISCUSSION 


We found the inequality n^> n c holds well for the structure Fig. 5. The same estimation 
can be done also for the ablative thermal protection structure 1 mentioned in the Introduc- 
tion. Note also that of the percolation threshold Vf c = ird 2 ln c /4:, which has been estimated 
in Ref. 12 for l/d = 80, yields Vf c = 0.03, while typical volume concentration of the fibers in 
such structures is Vf = 0. 1-0.2. This means that Eq. (3) can be used for the TC estimation, 
and the result is in agreement with the experiments. 2,14 

More interesting is the effect of anisotropy. In Sec. II we found the simple estimation 
Eq. (2), which depends on the correlation lengths Dy and D±, but is independent on the 
connectivity A. This result seems surprising, but it can be understood from the following 
example. Let us consider the thermal conductivity of a cubic mesh of fibers. It can be 
estimated as the thermal conductivity of the A x A x A cube. The number of disconnected 
fibers within the cube is (A/D) 2 and the resistivity of each fiber is, 4\/ird 2 Kf. The resistivity 
of the cube can then be written as (AX/nd 2 • (7tD 2 /4A 2 ) or 1/A K e ff- Thus, 



This result would not changed, if we had disconnected some of the contacts in the cubic 
mesh. The result changes only if we disconnect all the contacts along the fibers. Therefore, 
the error of the estimation Eq. (3) cannot be better than X/l ~ n c /n. This is the effect 
of dead ends, which is not important far from the percolation threshold n c . For the same 
reason, the small bridges of the matrix material do not affect the thermal conductivity in 
the case of the immersed fibers. 


Another reason for the error in Eq. (3) comes from the dispersion of resistivity of the 
fiber segments between two contacts. This effect can be estimated with the factor 21 

/ 2a 2 ' 


-l 


$ 


1 + 


N 


where a is the relative dispersion of the resistivity (i.e., relative dispersion of the length 
between adjacent contacts in our case) and N > 6 in the 3D resistor net. Assuming a < 0.5, 
we estimate the relative error as 0.1. 
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Typically, the matrix material (phenolic resin) can be porous itself. The size of its 
pores is very small, and their effect can be estimated, if we make the substitution 17 km — > 
km (l — | Pm) , where pm is the intrinsic matrix porosity. This makes the effect of the matrix 
inessential. 

Estimation of the conductivity of anisotropic structures can be performed using Eq. (3). 
Indeed, the fiber density Vf can be obtained from the structure porosity, while the ratio of 
the correlation lengths, D±/D\\, can be estimated by analyzing the oscillations of the Fourier 
transforms of the CF across and along the structure [Fig. 5(b)]. 

In order to calculate the correlation function, we have used the model structures. The 
same calculation can be made also using the real SEM or optical images of the cross-sections. 
Similar images can be obtained from the X-ray tomography [Fig. 5(a)], Basically one slice 
is enough for the CF calculation; however its accuracy can be significantly increased, if we 
use a few slices and average the results of calculations. 

In conclusion, we have analyzed the expressions (2) for the thermal conductivity of the 
fiber structures and investigated the estimation of their parameters from the information 
provided by the CF. 

This work was supported by the Swiss National Science Foundation (grant 200021 — 
130274/1). 
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